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Abstract 

This study describes the establishment of similarity conditions between two structural 
systems. Similarity conditions provide the relationship between a scale model and its 
prototype, and can be used to predict the behavior of the prototype by extrapolating 
the experimental data of the corresponding small scale model. Since satisfying all the 
similarity conditions simultaneously is difficult or even impossible, distorted models 
with partial similarity (with at least one similarity condition relaxed) are more practi- 
cal. Establishing similarity conditions, based on both dimensional analysis and direct 
use of governing equations, is discussed and the possibility of designing distorted mod- 
els is investigated. The method is demonstrated through analysis of the cylindrical 
bending of orthotropic laminated beam— plates subjected to transverse line loads. 
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Chapter 1 

INTRODUCTION 


1.1 Introduction 

Before production, any new design is subjected to many investigations through the- 
oretical analyses and experimental verification. As a system becomes more complex, 
assumptions are usually made in order to formulate a mathematical model for the 
system . In the absence of a complete design base, a new system requires extensive 
experimental evaluation until it gains the necessary reliability and desired perfor- 
mance. Since most of these tests are destructive, many test specimens are needed 
(Morton, 1988). 

For large and “oversize” systems, such as tall buildings , dams, bridges, spacecraft, 
airplanes, and space stations, creating the actual working conditions for testing the 
prototype most of the time is impossible, as in providing a “zero gravitional accelara- 
tion condition ” on the ground for testing large space stations or antennas (Letchworth 
et al., 1988 and McGowan et al,1990). Even when a prototype test is possible, it is 
expensive, time consuming , and difficult to control. Thus, it is extremely useful if 
a prototype can be replaced by a similar scale model which is much easier to work 
with. The only possible way to obtain experimental data of overall performance of 
such a system and the interaction of its elements is to design a small similar system 
( scale model ) which replicates the behavior of the actual system ( prototype ). The 
accuracy of the behavior of the prototype , which is predicted from interpreting the 
test results of the model, is dependent on the relationship between the corresponding 
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variables and parameters of model and its prototype. 

Similarity of systems requires that the relevant system parameters are identical 
and these systems are governed by unique set of characteristic equations. Thus, if 
a relation or equation of variables is written for a system, it is valid for all systems 
which are similar to it (Kline, 1965). Each variable in a model is proportional to 
the corresponding variable of the prototype. This ratio, which has essential roles in 
predicting the relationship between model and its prototype, is called scale factor . 

Since making a precise and adequate experiment is difficult and expensive, it is 
more convenient to run a series of experiments in nondimensional form of variables, 
the results of which can then be used for a similar system. Models, as a design aid, 
have been used for several years, but the use of scientific models which are based on 
dimensional analysis was first discussed in a paper by Rayleigh. Similarity conditions 
based on dimensional analysis have been used since Rayleigh’s time (Macagno,1971), 
but the applicability of the theory of similitude to structural systems was first dis- 
cussed by Goodier and Thomson (1944). They presented a systematic procedure for 
establishing similarity conditions based on dimensional analysis. 

In establishing similarity conditions between the model and prototype two proce- 
dures can be used. The similarity conditions can be established either directly from 
the field equations of the system or, if it is a new phenomenon and the mathematical 
model of the system is not available, through dimensional analysis. In the second 
case, all of the variables and parameters which affect the behavior of the system must 
be known. By using dimensional analysis, an incomplete form of the characteristic 
equation of the system can be formulated. This equation is in terms of nondimen- 
sional products of variables and parameters of the system. Then, similarity conditions 
can be established on the basis of this equation. 

The objective of this study is to discuss and demonstrate the establishment of 
the similarity relations between two structural systems. Similarity conditions provide 
the relationship between model and its prototype, and can be used to extrapolate 
the experimental data of a small and less expensive model in order to predict the 
behavior of the prototype. In all of our work in this area we will restict ourselves to 
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linearly elastic material behavior. 


1.2 Literature Review 

The scientific small scale model based on dimensional analysis, was first discused by 
Rayleigh(1915). He established the fundamentals of dimensional analysis based on 
Fourier’s work. This principle has been reviewed and completed by Riabouchinsky, 
Buckingham, Bridgman, Bickhoff, and Langhaar. (for more detail see Macango, 
1971). 

The applicability of the theory of similitude to structural systems was first dis- 
cused by Goodier and Thomson(1944) and later by Goodier in 1955. In the 50’s 
and 60’s many interesting books have been published in this area, Murphy(1950), 
Langhaar(1951), Sedov(1959), Kline(1965), Skoglund(1967). Most of these authors 
discussed similitude theory based on dimensional analysis. Kline(1965) gives a per- 
spective of the method based on both, dimensional analysis and the direct use of 
the governing equations. Szucs (1980) is particularly thorough on the topic of the 
similitude theory. He explains the method with emphasis on the direct use of the 
governing equations of the system. 

Many research activities have been conducted on modeling of dynamic and static 
behavior of structural systems, especially on modeling the reinforced concrete struc- 
tures, (see Sabines and White (1966,1977), Harris et.al.(1966,1970)). Krawinkler 
et.al.(1978) described detailed model studies on earthquake resistance of structures 
and presented a comprehensive reference on the dynamic behavior of structures for 
seismic engineering analysis. 

Since reinforced composite components require extensive experimental evaluation, 
there is a growing interest in small scale model testing. Morton (1988) discusses the 
application of scaling laws for impact-loaded Carbon-Fiber composite beams. His 
work is based on dimensional analysis. Qian et.al.(1990) conducted experimental 
studies of impact loaded composite plates, where the similarity conditions were ob- 
tained by considering the governing equations of the system. 
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In recent years, due to large dimensions and unique structural design of the pro- 
posed space station, small scale model testing and similitude analysis have been con- 
sidered as the only option in order to gain experimental data. Letch worth et.al.(1988), 
Shih et.al.(1987), Hsu et.al.(1989), and McGowan et.al.(1990) discussed the possibil- 
ity of scale model testing of space station geometries especially for vibration analysis. 

Most of these studies have used complete similarity between model and prototype. 
The objectives of the investigation described herein are: 

• explore two fundamental methods of similitude analysis 

• create necessary similarity conditions in order to design an accurate distorted 
model 

• evaluate the derived similarity conditions analytically and correlate the actual 
experimental data of the prototype with the small scale model predictions. 

The experimental data of cylindrical bending of several beam-plates are used to sup- 
port the applicability of the derived similarity conditions. 
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Chapter 2 

THEORY OF SIMILITUDE 


This chapter is devoted to consideration of foundations of the similitude theory. The 
meaning of similitude is explained, and two major methods of obtaining similarity 
conditions are discussed. 

Similitude theory is concerned with establishing necessary and sufficient conditions 
of similarity between two phenomena. Establishing similarity between systems helps 
to predict the behavior of a system from the results of investigating other systems 
which have already been investigated or can be investigated more easily than the 
original system. Similitude among systems means similarity in behavior in some 
specific aspects. In other words, knowing how a given system responds to a specific 
input, the response of all similar systems to similar input can be predicted. 

Euclid established the fundamental theory of similitude by defining geometric sim- 
ilarity for plane geometric figures. According to this theory, enlarging or contracting 
all dimensions of a figure by a constant ratio, which is called scale or similarity factor, 
forms figures which are similar to the original figure. 

The behavior of a physical system depends on many parameters, i.e. geometry , 
material behavior, dynamic response, and energy characteristic of the system. The 
nature of any system can be modeled mathematically in term of its variables and 
parameters. A prototype and its scale model are two different systems with different 
parameters. The necessary and sufficient conditions of similitude between prototype 
and its scale model require that the mathematical model of the scale model can be 
transformed to that of the prototype by a bi-unique mapping or vice versa (Szucs, 
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1980). It means, if vectors X p and X m are the characteristic vectors of the prototype 
and model , then we can find a transformation matrix A such that 


X p = \X m or X m = \~ l X p (2.1) 

The elements of vector X are all the parameters and variables of the system. A 
diagonal form of the transformation matrix A is the simplest form of transformation. 
The diagonal elements of the matrix are the scale factors of the pertinent element of 
the characteristic vector X . 


A x i 0 . . . 0 

0 A X 2 ... 0 


0 0 ... A xn 


( 2 . 2 ) 


where X X) — — denotes the scale factor of a:,-. In general the transformation 

^im 

matrix is not diagonal. 

Since similitude theory gives many alternative ways for investigating a system, it 
has been used in areas which primarily involved many experimental investigations. 
Two methods, dimensional analysis and direct use of governing equations are used to 
establish similitude between systems. 


2.1 Dimensional Analysis 

The principal purpose of dimensional analysis is to reduce the number of parameters 
by establishing product groups of variables (t t- terms) such that all terms are mutually 
independent and dimensionless. 

Rayleigh established the fundamentals of dimensional analysis based on Fourier’s 
work. Riabouchinsky proved the fundamental theorems of dimensional analysis , 
and Buckingham reformulated Riabouchinsky’s theorems. Buckingham’s theorem 
was discussed and revised by Bridgman, BrickhofF, Langhaar, Van Driest, and Brand, 
and was named Buckingham’s 7r-Theorem, or simply the 7r-Theorem (Macagno,1971). 
This theorem states that a given relation among n variables can be reduced to an 
equivalent relationship in terms of a complete set of n — r dimensionless 7r-terms, 
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where r is the rank of the dimensional matrix. 


<f>(x u x 2 ,x 3 , •••,!„) = 0 or $(7Ti,7r 2 , •••,7r n _ r ) = 0 (2.3) 


The function $ is called functional relation . The functional relation need not be a 
known function of the system variables. Function $ only shows the reduced form of 
relevant variables and does not give any information about the nature of the solution. 
Any physical system can be described in term of various combinations of fundamental 
quantities such as M(mass) , L(lenght) , F(force) , T (time), and secondary quantities 
such as v (velocity) ,<x (stress), A(area) etc. The first step is to choose an adequate 
set of fundamental quantities (dimensions). Then all variables can be described in 
this system, which is unique for that set of fundamental dimensions, and can change 
by choosing another set of dimensions. It is necessary that all pertinent parameters, 
even constants like the gravitational acceleration be included, otherwise the analysis 
leads to incorrect and incomplete results. 

Barr(1983) has presented five methods for dimensional analysis. These are the 
Rayleigh method and various modifications such as the Buckingham method, the 
Echlon Matrix method, the Basic Stepwise method, and the Proportionalities method. 

As an example consider the procedure for finding the buckling load of a simply 
supported elastic column with an elastic support at the middle (Figure 2.1). The 
relevant variables are: 

variable dimensions 


w=deflection L 

L=span L 

E=modulus of elasticity F L ~ 2 
I=moment of inertia L 4 

k=$pring stiffness F L~ l 

P=axial load F 





Figure 2.1: Simply Supported Column 


7 



<f>(w, L, E , /, k, P) = 0 

The dimensional matrix and the resulting 7r-terms are [see Barr (1979) for details on 
the procedure]: 

L E w I k P L E w I k P 

L ( 1 -2 1 4 -1 0\ l/l0 1412\ 

F^O 1 00 1 1 J ^ 1 0 0 1 1 j 

w I k P 

7fl “ L ,7r2_ L 4 ,7F3_ LE ’ ** ~ EL 2 

By applying one of the dimensional analysis methods, a set of tt- terms is calculated 

which is not unique and many other combinations of parameters can yield several sets 

of 7r-terms. 

Since experimental control of some variables is easier than others, it is more con- 
venient to have these variables only in one of the 7r-terms. Power products of 7r-terms 
are also nondimensional which can be used in place of the inadequate tt- terms. The 
resulting 7r-terms can be written as 

^(*1, 7T 2 , 7T 3 , 7T 4 ) = 0 

If 7Tj involves the dependent variable, it can be written as a function of the other 
7r-terms (independent variables). 

TT! = $(7r 2 , 7r 3 ,7r 4 ) 

Similarity conditions require that the equations of two similar systems be the 
same. If the 7r-terms of the functional equation for two systems are the same, then 
= $2 even if we don’t know the functional equation completely. These equalities 
of 7r-terms which determine the conditions for which the two systems are similar, 
are called the similarity conditions or scaling laws for these systems and for specific 
phenomena. 

^1 p ^p(^2p> ^3p> ^4p) i ^lm (^2m) ^~4m ) 

if TTtm = n ip f or * = 2, 3, 4 , then 

^p(^*2p) ^*3p> ^*4p) = ^m(^2m)^3mi ^4m) 
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hence 


^"lm — ^*1 p 

Since these 7r-terms are combinations of geometric , dynamic, material, and kine- 
matic parameters of the systems, the above equalities define different similarities, 
such as geometric , material, kinematic, and dynamic similarity (Schuring,1977). 

2.2 Direct Use of Governing Equations 

The field equations of a system with proper boundary and initial conditions char- 
acterize the behavior of the system in terms of its variables and parameters. The 
dependent variables or response of the system is a function of its independent vari- 
ables (input and parameters of the system). If the field equations of the scale model 
and its prototype are invariant under transformation A and A -1 , then two systems are 
similar [see Eq.( 2.1)]. This transformation defines the scaling laws (similarity condi- 
tions) between all parameters, input, and response of the two systems. The problem 
of finding similarity conditions for the buckling of a column, which was analysed 
by dimensional analysis, is considered again (Figure 2.1). Considering the symmet- 
ric mode of deflection, the vertical reaction of spring at the middle is k8 . Where 
8 is deflection at the middle of column. From symmetry, for half of the column 
0 < x < L/2 the reaction force is k8/2 . From Simitses (1976) the differential 
equation and associated boundary conditions of the system are: 


rvu u/ 

EIw <xx + Pw = -y 

' u;(0) = 0 

B.C.l w tX (L/2) — 0 
, w{L/2) = 8 

For model and prototype we may write 

F i , P „, _ K to 

EpIp dxj + Pp p ~ 2 


E m I m 

m m dx 2 


j-. k m 8 m W m 

— + PmWrn = 

m " 


where subscripts m and p refer to model and prototype respectively. 
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By defining scale factors A i , the variables of the prototype can be written as 
Zip = A Xi Xi m . The similarity conditions between model and prototype are determined 
by substitution of the A x ,Xi m into the differential equation of the prototype and by 
requiring that the result be the differential equation of the model, Eq.( 2.5)(complete 
similarity). 


^AtfA/Au,^ p J <PWrn , ( ^ ^ \T3 ... / \ \ \ 

V ^2 / * m ^^,2 ' \^p^w)*rn^m — 2 


(2.7) 


Eqs.( 2.6) and ( 2.7) are the same if the terms in parentheses of Eq.( 2.7) are all equal. 

(2.8) 


(^^) = (A„A„) = (A t A s A„) 


AJ 

Dividing Eq.( 2.8) by one of these terms, i.e. the first term , 

J A,a; AtAtAjS 

Xe^i X e Xi 


or 


and 


a p a; 


A^A/ 


= l 


= l 


or 


A^A/ 

Similarly for the boundary conditions, 

Am _ 

a7“ 


or 


A p = \ x 2 \ e^i 




or 


A,, r — At 


(2.9) 


(2.10) 


( 2 . 11 ) 


( 2 . 12 ) 


Eqs( 2.10)— ( 2.12) are necessary and suficient conditions for complete similarity 
between model and prototype. These similarity conditions define three relation- 
ship among seven unknown (A’s) . Hence we can choose four of A’s freely and 
find the values of the others by requiring satisfaction of the similarity conditions, 
Eqs.( 2.10) — ( 2.12). 
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Chapter 3 


PARTIAL SIMILARITY AND 
DISTORTED MODEL 


If all the 7r-terms for the model and the prototype are the same, or if all similarity 
conditions are satisfied, the two systems are completely similar. But often complete 
similarity is difficult or even undesirable. The model which has some relaxations in 
similarity conditions is called a distorted model. 

In complete similarity (in terms of 7r-terms) tt i m = 7 r; p for i = 1, • • • , N but 
in partial similarity 7 r tm = 7r ip for i = 1 , * • * , k where k < N . These relaxations 
in the relationship between two systems cause model behavior to be different from 
that of the prototype. Understanding these relaxations and their effects in model 
behavior can be used to modify the model test data so as to predict the behavior of 
the prototype. 

3.1 Dimensional Analysis 

In dimensional analysis, similarity of each 7r-term of the model with corresponding 
7 r-term of the prototype is a function of the scale factors of the variables. Since there 
are m — r 7r-terms, the results will be m — r functions of m scale factors. If the 
two systems are completely similar, the number of scale factors which can be chosen 
freely is equal the rank of the dimensional matrix, r , and the m — r remaining scale 
factors are determined by solving m — r similarity conditions. 
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complete — similarity 


partial — similarity 


^lm — p 


^lm — ^*1 p 


^*2m ^2p 


^2m ^*2p 


—7r lp ^fcp 

7r fc+lm = 7r /:+lp TTAr+lm ^ ^Ar+lp 

7r nm =7r np ^nm 7^ ^np 

For a distorted model, inequality of some 7r-terms causes a change in scale factors 
of those terms. 7r-terms which include these scale factors will also change. In this case, 
the number of unknown scale factors is greater than the number of similarity equa- 
tions and additional relationships between variables are needed. Dimensional analysis 
cannot provide these relationships. These additional relationships between variables 
can be established by the governing equations of the system, such as equations of 
equilibrium and compatibility, kinematic relations, material behavior equations and 
boundary conditions. If these equations are not available, by conducting a series of 
experiments the effect of distortions on each term can be found while the other terms 
are kept constant. A sufficient amount of data should be determined from model tests 
so that these relationships can be understood clearly (Langhaar,1954). 


3.2 Direct Use of Governing Equations 

When governing equations of the system are used for establishing similarity condi- 
tions, the relationships between variables are forced by these equations. Suppose 
the system has m variables and similitude analysis of the governing equations of the 
system define n relationship among m unknown, (scale factors of these variables). 
If the two systems are completely similar m — n scale factors can be chosen freely 
and the values of the other scale factors are found by using n similarity conditions. 
When at least one of the similarity conditions can not be satisfied, partial similarity is 
achieved. In this case, since each variable has different influence on the response of the 
system, the resulting similarity conditions have different influence. By understanding 
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the effect of variables and similarity conditions over desired intervals, the similar- 
ity conditions which have the least influence can be neglected without introducing 
significant error(Kline,1965). 

Suppose we want to design a reasonable (able to test) model for a large rectangular 
plate. The plate is simply supported at all edges and loaded with uniform load 
of intensity q. Assume uniform cross section and isotropic material, the governing 
differential equations and boundary conditions are well known (Timoshenko, 1959) 

d 4 w , 2 d 4 w d 4 w = q 
dx 4 dx 2 dy 2 dy A D 


and B.C at x = 0, a 


w = 0 

Mx = ~ D ltf = ° 


and at y — 0, b 


w = 0 
d?w 

M, = -Z>-=0 


By applying similitude theory 

^<7 / ry r \ 

a|a* = K ~ ^ 

Now to find the scaling laws from Eq.( 3.5), we have three choices. Dividing Eq.( 3.5) 
by first term, yields 

A, = A, , A. = ^ (3.6) 


Dividing Eq.( 3.5) by the second term, yields 


A x = A 


y > ''w 


A,Ap; 
— % 


and finally dividing Eq.( 3.5) by third term 


Ajp — A 


y ) 


. _A,AJ 

\ 
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To find which one of Eqs. ( 3.6)-( 3.8) gives the best prediction for the prototype 
behavior, the theoretical deflections of the model are projected with each condition 
and compared to theoretical deflection of the prototype. For complete similarity, all 
three give the same results. In order to avoid impractical size for the cross section we 
need to choose different scale factors in the x ,j/, and z directions. 

Suppose A y = C\ x where C > 0. Eqs.( 3.6)-( 3.8) can be written as 


\ _ _ /V Wf 

A W X ^ X 

a,a;a; A,Aj 

X ^ X 

. Vj pA,Aj 

S\ w Kj . 


c = 1 
c = c 2 
c = c 4 


(3.9) 

(3.10) 

(3.11) 


To find which one of Eqs.( 3.9) — ( 3.11) gives the best prediction for the prototype 
behavior, the theoretical deflections of the model are projected with each one of 
these conditions, Eqs.( 3.9) — ( 3.11), and compared to the theoretical deflection of 
the prototype. 

There are three possible configurations for the prototype. 

• Rectangular 


a 

b >l 


• Square 

• Inverse rectangular 


a 

1 = ' 


a 

6 < 1 


The characteristics of the model were calculated by using the given scale fac- 
tors. The theoretical deflections of the model and prototype are given by the same 
expression. From Timoshenko (1959), the deflection is given by 


16<7o v— ^ 


. rmrx . niry 
oo oo sin szn- 


Jl 


t t- 6 7~) ^ r/ ^ \2 / ^ \2i2 

* U rn=l n=l mn [( — ) 2 + ( y ) 2 ] 2 
a b 


(3.12) 


By using Eqs.( 3.9)— ( 3.11) the experimental data of the model were projected to 
predict the deflections of the prototype. The predicted deflections of the prototype 
were compared to the theoretical deflections of the prototype. The result is the % 
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discrepancy between predicted and theoretical deflections of the prototype. 


%Discr. = 


\wth. - *y.| 

Wth. 


x 100 


(3.13) 


The % discrepancies of the prototype are calculated for different values of C (dif- 
ferent A y ). The different values of C imply different aspect ratios (j) for the model. 

b 

For each C the experimental data were projected by Eqs.( 3.9),( 3.10), ( 3.11). The 
possibility of distortion in the y-direction for each case is investigated. In this analy- 
sis the experimental data of the model are manufactured from theoretical deflections 
of the model by randomly introducing a ±10% discrepancy. For all of these cases 
A e K A« 1. 


Case— 1 : 

When the prototype is a rectangular plate (-r)p > 1, a wide range of models can be 

a o 

used. In this case as (7 — ♦ 10 , (— ) m increases, Eq.( 3.11) yields excellent predictions 
and the % discrepancy for the prototype decreases. Similarly, when C —* 0.01 the % 
discrepancy increases. In this case, Eqs.( 3.9),( 3.10) are not good choices for similarity 
condition. In other words a rectangular plate can be replaced by another rectangular 
plate with different aspect ratio, where (t)p — 1 a nd Eq.( 3.11) provides the needed 
similarity condition. Figures 3.2— 3.7 present the predicted and theoretical deflection 
of the rectangular plates for different values of C. 

Case— 2: 


In this case the prototype is a square plate. The distortion is restricted to values 
of C which are close to one. It means that the model must be a rectangular with 
0.83 < ( — ) m < 1.2. For this case Eqs.( 3.9),( 3.10) yield good result. As C — > 10 or 
C — * 0.01 the %discrepancy increases rapidly and none of the similitude conditions, 
Eqs.( 3.9),( 3.10),and( 3.11), yield acceptable results. In this case Eq.( 3.11) is not a 


,a. 


suitable similarity condition even for small distortion((-) m « 1). Figures 3.8— 3.10 

b 
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show the predicted deflection of the square plates for different values of C. 

Case— 3 : 

In this case where the prototype is a rectangle with (^)m < 1 (Inverse rectangular), 
Eqs.( 3.9) predicted the behavior of prototype very well, when C — * 0.01. However 
as C —* 10 the %discrepancy increases slowly. The similarity conditions Eqs.( 3.10) 
and ( 3.11) are not suitable at all, since distortion from complete similarity causes the 
%discrepancy to increase very fast. Figures 3.11— 3.15 demonstrate the predicted 
deflection of the inverse rectangular plates for different values of C. 
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Table 3.1: Characteristics of the Prototypes 


prototype 

a (in) 

b (in) 

h (in) 

E (10 6 psi) 

V 

Rectangular 

6000 

600 

60 

25 

0.23 

Square 

6000 

6000 

60 

25 

0.23 

Inverse rect. 

3000 

6000 

60 

25 

0.23 


Table 3.2: Characteristics of the Models 
m °d e l C a (in) b (in) a/b h (in) 


R1 

0.01 

30 

R2 

0.05 

30 

R3 

0.1 

30 

R4 

0.5 

30 

R5 

1 

30 

R6 

1.5 

30 

SI 

0.1 

30 

S2 

0.5 

30 

S3 

1 

30 

S4 

1.1 

30 

S5 

1.5 

30 

S6 

2 

30 

S7 

10 

30 

IR1 

0.1 

15 

IR2 

0.5 

15 

IR3 

0.75 

15 

IR4 

1 

15 

IR5 

2 

15 

IR6 

4 

15 

IR7 

10 

15 


R : Rectangular 

S : Square 

IR : Inverse rectangular 


2 

300 

0.1 

0.3 

10 

60 

0.5 

0.3 

20 

30 

1 

0.3 

100 

6 

5 

0.3 

200 

3 

10 

0.3 

300 

2 

15 

0.3 

20 

300 

0.1 

0.3 

100 

60 

0.5 

0.3 

200 

30 

1 

0.3 

220 

27.3 

1.11 

0.3 

300 

20 

1.5 

0.3 

400 

15 

2 

0.3 

2000 

3 

10 

0.3 

20 

300 

0.05 

0.3 

100 

60 

0.25 

0.3 

150 

40 

0.625 

0.3 

200 

30 

0.5 

0.3 

400 

15 

1 

0.3 

800 

7.5 

2 

0.3 

2000 

3 

5 

0.3 
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Figures 3.16— 3.18 present the discrepancies between theoretical and predicted 
maximum deflections of the simply supported plates as a function of C, the cross 
section coefficient of the model, where 0.001 < C < 10. Each C represents a specific 
A y or \ z . 

Figure 3.16 shows these discrepancies for rectangular plates when different sim- 
ilarity conditions are used. It is observed that Eq.( 3.11) can predict the maximum 
deflection of the model very well for a large range of C; especially for C > 1. The 
other conditions, Eqs.( 3.9)and ( 3.10), have a large discrepancy for small distortion 
of the model. 

Figure 3.17 presents the discrepancies for square plates. It is clear that none 
of Eqs.( 3.9)— ( 3.11) are suitable representation of similitude for distorted models. 
However, for small changes in C (i.e where C is close to 1) Eq.( 3.10) is reasonably 
accurate. 

The discrepancies for inverse rectangular plates are plotted in Figure 3.18. In 
this case Eq.( 3.9) can predict the behavior of the prototype very well as long as the 
configurations of model and prototype are the same. 
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Figure 3.2: Predicted , Eq.(3.11), and Theoretical Deflection of S-S Plate , (— ) p > 1 
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Figure 3.8: Predicted, Eq.(3.10), and Theoretical Deflection of S-S Plate , (— ) 
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Figure 3.10: Predicted, Eq.(3.11), and Theoretical Deflection of S-S Plate, ( 
for C = 0.1 . 
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Chapter 4 

APPLICATION TO 
CYLINDRICAL BENDING OF 
LAMINATED BEAM-PLATES 


In this chapter, as an initial effort, similarity conditions are developed in order to 
design resonable, distorted scale models for orthotropic laminated beam-plates. Plates 
are subjected to transverse line loads. Later, the experimental data from 3— point 
tests ( cylindrical bending) of the plates are used to verify the accuracy of the distorted 
model in order to predict the maximum deflection of the prototype. 

The available experimental data is used in the following way. One of the plates is 
considered as prototype, the other as its scale model. Then the data of the model are 
projected by scaling laws in order to predict the prototype behavior. The predicted 
data are compared to actual experimental data of the prototype. 

4.1 Similitude Analysis Based on Dimensional 
Analysis 

Consider a cross-ply laminated plate composed of N orthotropic layers. The plate is 
simply supported at x = 0, a, free at y — 0, b ,and subjected to a transverse line load 
with intensity q 0 (Fig. 4.1). The deflection of the system can be written as function 
of these variables ; 
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10 = deflection 

L 

a = span 

L 

Eu — long. mod. of elas. 

FL~ 2 

E22 — trans. mod. of elas. 

FL~ 2 

G 12 — modulus of rigidity 

FL~ 2 

9 =. fiber orientation 

none 

N — number of lamina 

none 

h = thickness 

L 

b = width 

L 

U12 = poisson’s ratio 

none 

^21 — poisson’s ratio 

none 

q Q = line load 

FL- 1 


<f>{w, a , En i £22, Gn, V\ 2 , V21 1 h, b , A 1- , <?) — 0 
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The result 7r-terms are 


w a b 

n i — — = T i ^3 — r 

a h h 

*6 = i/ 2 i ,7 t 7 = N ,7 r 8 = 0 

where functional relation can be written as 


, 7T 4 — 


, 7Tg — 


£„ 

E 22 

E „ 


,7T 5 — J/ 12 


i^IO — 


90 


£11 ft 


12 


w _ a b G 
a = ’ft 


at /) E 22 qo v 

,*n ,N ,0 ,— ) 


complete similarity is achived if 


^l<21 ^0 1 

^J?ll ^£22 ^(*12 

Au- = A(, = A/, = Ajr, 

Assume the model and prototype have the same material and fiber orientation. Then 
from dimensional analysis (complete similarity) 

^12 ^*^21 ^i?ll ^&22 ^Gl2 ^0 — 1 (4.1 ) 

^ w — = ^ h (4.2) 

If different scale factors are used in the x,y,z directions. Then A a ^ X b / \ h 
hence A,*, ^ A a and a new relationship for A„, must be found. Dimensional anal- 
ysis cannot provide any additional information. This necessary information can be 
extracted from the analysis of laminated plates. Since model and its prototype both 
have identical governing differential equations, similarity of these equations gives the 
additional relationship. Thus, we have to apply similutude theory through the gov- 
erning equations. For this reason in partial similarity, it is more convenient to employ 
similitude theory based on direct use of field equations in order to design a reasonable 
and acceptable distorted model. 


30 



4.2 Similitude Analysis Based on Direct Use of 
Governing Equations 


By assuming that the displacement functions are independent of y , or u = u(x) , 
v = 0 , w = w(x) (cylindrical bending) , from Ashton and Whitney (1970) the 
governing differential equations and boundary conditions are reduced to 

d 4 w qAn 

~d- * ~ A n D u - Wi (4 ‘ 3) 
cPu B n d 4 w 

lx 3 = Aldl* ( 4-4 ^ 

and the B.C.’s at x — 0, a are 


w = 0 

(4.5) 

. du _ d?W 

N X x — An j Bn _ 0 

ax ax 1 

( 4 . 6 ) 

, , du cPw 

M XX B\\ j D\\ — 0 

ax dx l 

(4.7) 

Eq.( 4.3) can be written as 


(A n D„ - = qA n 

( 4 . 8 ) 

By applying similitude theory, the resulting similarity conditions are 


Aj4h ^Dn A^Ag 

( 4 . 9 ) 

or 


^A n A Du = A| n 

(4.10) 

A^A A^A^ 

(4.11) 

Similarly from Eqs.( 4.4),( 4.6), and ( 4.7) we have 



(4.12) 

Abu A U A X — A^A^jj 

( 4 . 13 ) 


The condition depicted by Eq.( 4.13) does not give any new information, since it can 
be obtained by combining Eqs.( 4.10) and ( 4.12). So, Eqs.( 4.10) — ( 4.12) denote the 
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necessary conditions for complete similarity between scale model and its prototype. 
For better understanding the restrictions of Eq.( 4.10), consider the definition of 
and D mn . 

N 

Amn ~ y!(Qmn)j(^7 — z j- 1) 

B m » = \ E(<5„n),(^ - *’_,) 

z i=i 
^ i=i 

where Zj is the coordinate of the upper surface of the j th lamina (measured from the 
plate reference surface). Let zj = Cjh where —0.5 < Cj < 0.5 and h is the total 
thickness (j = 0, 1, • • • ,N). (Q mn )j, the transformed stiffnesses for the j th lamina is 
given in terms of the engineering orthotropic constants and the fiber orientation angle 

e. 

Q — f( 0 , E\u E22, ^ 12 , G12) 

This allows us to express A mn ,B mn , D mn in terms of h and functions of all <5,./V,and 
the sequence in which the plies are arranged. 



A m n — hfa(Qmm -^0 



Bmn = h 2 f b (Q mn ,N) 
D mn = h 3 f d (Q mn ,N) 


or as scale factors 

^A mn = ^hFa(Qmn,N) 

(4.14) 


*B mn =\lF b (Q mn,N) 

(4.15) 


*D mn = \lF d (Q mn ,N) 

(4.16) 

where f = m»k 

fi{Q,N) m 

Substituting Eqs.( 4.14) 

yl Qj by d. 

— ( 4.16) into Eq.( 4.10) we have 



F a {Qu, N)F d (Q U y N ) = F 2 (Q U , N) 

(4.17) 
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This condition, Eq.( 4.17), is satisfied if the model and prototype are made of the 
same material with identical N and the same stacking sequence of the laminae. 


4.3 Analytical Verification 

In this section, the accuracy of the derived behavioral similarity conditions, Eqs.( 4.10) 
— ( 4.12) is evaluated analytically, in order to determine the level of confidence that 
can be expected in interpreting the data from the model experiments. 

Consider a cross-ply laminated E-Glass/Epoxy plate composed of 96 orthotropic 
layers (0/90/0/.. .) 96 as the prototype. We desire to find the maximum deflection 
of the prototype by extrapolating the pertinent values of a small scale model. The 
model has the same stacking sequence as the prototype but with a smaller number of 
layers. The prototype and its scale model have the following characteristics; 


prototype (0/90/0 .. .) 96 : 

a — 90 in 

b — 100 in 

h = 0.858 in 

N =96 

model (0/90/0 . . ,) 16 : 

a = 5.0 in 

b = 6.139 in 

h = 0.143 in 

N = 16 

Scale Factors : 

X a = 18 

X b = 16.29 

II 

S© 

II 


For simplification we assume that model and prototype have the same material 
properties (A.e u = Xe 22 = A„ 12 = 1), and A, = A*. By employing similarity condition 
Eq.( 4.11), (note that \ P = A^A,; therefore the condition becomes X w Xd u = A^A p) 
the theoretical maximum deflections of the model are projected in order to predict 
the maximum deflections of the prototype. Figure 4.2 presents the theoretical and 
predicted maximum deflections of the prototype and corresponding theoretical values 
of the scale model. The derived scaling laws can be used with high level of accuracy in 
predicting the prototype behavior. Note that the model was designed by employing 
the free scaling factors. 
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Locd (lb/in) 


Figure 4.2: Theoretical and Predicted Maximum Deflections of Prototype (0/90/0/... )s$ 
when model (0/90/0/. ..)i 6 is used, (\g u = A e 22 = A„ 12 = 1, A a = 18 , X b = \ q - 
16.92 , X h = \ N = 6). 


If the model and prototype have the same material properties then Ag = 1 • For 
a scale model which has equal number of plies and stacking sequence of laminae 
with prototype F,(Q mn ,N ) = 1, i = a,b,d, and therefore A j4lJ , A 5ll , and X Dll are 
equal to h , /i 2 ,and h 3 respectively. In this case the similarity condition, Eq.( 4.10), is 
automaticly satisfied and Eqs.( 4. 1 1),( 4.12) can be written as 

A w = A 2 A p A a ~ 3 (4.18) 

X u = X w X h X x 1 (4.19) 

In general, by choosing the model material and using Eq.( 4.10), the number of 
plies of model ( N m ) can be determined. Since N m must be an integer number, it is 
difficult to satisfy the Eq.( 4.10) , therefore partial similarity with a distorted model 
is achieved. 

Consider that model and prototype have the same material properties with dif- 
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ferent number of layers N m ^ N p . Having material properties for both systems and 
N p , Eq.( 4.10) determines an approximate value for N m . Using these parameters 


where Cj 


with c a = 


, KMq,n) p 

D " hlMQ,N) m ° dXh 

(4.20) 

- f d (Q,N) r „ “d EqS '( 41°)-( 412 > >™ ld 


= Cj 1 A x 3 A/ l _3 A/> 

(4.21) 

A u A w A x Ah 

c a 

(4.22) 

A u A w A x Ah 

Cb 

(4.23) 

MQ,N) r MQ,N )„ 

UQ,N) m ’ ‘ MQ,N) m 
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4.4 Experimental Verification 


To demonstrate the use of the above analysis, the experimental data of cylindri- 
cal bending tests of 10 orthotropic plates are considered. These plates are made 
of E-Glass/Epoxy and Kevlar/Epoxy with different number of layers and stacking 
sequences of the laminae. Experimental data were provided by Professors Shive 
Chaturverdi and Robert Sierakowski (1991). 


Table 4.1: Characteristics of the Employed Beam— Plates 


Plate 

b(in) 

h(in) 

N 

sequence 

G x 

6.139 

0.143 

16 

(0/90/0... ) 16 

g 2 

5.975 

0.147 

16 

(0/90/0... ) 16 

G3 

6.087 

0.137 

15 

(O3/9O3/O3/9O3/O3) 

g 4 

6.101 

0.131 

15 

(O3/9O3/O3/9O3/O3) 

k 5 

6.109 

0.132 

18 

(0/90/0... )i 8 

K 6 

6.111 

0.147 

18 

(0/90/0... ) 18 

k 7 

6.042 

0.155 

20 

(O4/9O4/O4/9O4/O4) 

K 8 

6.108 

0.147 

20 

(O4/9O4/O4/9O4/O4) 

k 9 

6.033 

0.135 

18 

(0 6 /90 6 /0 6 ) 

K 10 

6.107 

0.142 

18 

(06/906/06) 


G : E-Glass/Epoxy 
K : Kevlar/Epoxy 


Table 4.1 shows the characteristics of the plates. All plates have identical span 
a = 5.0 in . In the demonstration one of these plates is considered as prototype and 
another one as its scale model. The derived similarity conditions, Eqs.( 4.18) — ( 4.23), 
are used to project the experimental data of the model in order to predict the deflec- 
tion of the prototype. By comparing the value of predicted deflection with the actual 
experimental deflection, the amount of discrepancy is calculated. The percentage of 
discrepancy is defined as 

%Discr. = ■ U>egp ' — x 100 (4.24) 

Wexp. 
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CASE - 1 


In this case, plate G4 is considered as the prototype and G3 is its scale model. 
From the data of Table 4.1 these scale factors are calculated as; 

A a = 1.0 , X b = 0.9977 , X h = 1.0458 

Since the number of the plies and the stacking sequence of the laminates are identical, 
then X w is [Eq.( 4.18)] 

A*, = X a 3 X q Xh~ 3 (4.25) 

Since A 0 and X b are known, Eq.( 4.25) relates A„, and A,. By choosing A, , the 
corresponding A^ is determined. A, must be chosen based on elastic limit of the 
model and the prototype. In this case, since both the model and prototype are made 
of the same material, A, can be expressed as a function of size scales (A a , A;,, or X b ) 
and A£ n . For different nondimensionalized loads, A, can be written as 


\q — j A{j “ 1 

(4.26) 

Ag A^ j A^ "■ A^ 

(4.27) 

Ag Xeh ^ h A fi 

(4.28) 


Figures 4.3 and 4.4 show the theoretical, predicted and experimental values of the 
maximum deflection of the prototype G4 and its model G3 as functions of load. As 
shown the model can predict the maximum deflections of the prototype with very high 
accuracy. For clearer presentation of the results, the % discrepancy between the model 
(experimental and theory) and its prototype (predicted and experimental , predicted 
and theoretical) are shown in Fig 4.5 and 4.6. The % discrepancy between predicted 
and experimental deflections of the prototype are less than 6%. This indicates an 
excellent prediction. The predicted data match very well with experimental and 
theoretical deflections of the prototype. 
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Load (lbs) 


Figure 4.3: Predicted and Actual Test Results of E-Glass/Epoxy Plate G3 

(O3/9O3/O3/9O3/O3) and its Model G{ (O3/9O3/O3/9O3/O3). 







Figure 4.6: %Discrepancy for Prototype G3 (O3/9O3/O3/9O3/O3) When Model Gf 

(O3/9O3/O3/9O3/O3) is Used. 
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CASE - 2 


In this case, two Glass/ Epoxy Plates G2 and G4 are considered as the prototype 
and its scale model respectively. From the data of Table 4.1, the scale factors are 
calculated as; 


A a = 1.0 , A 6 = 0.9793 , X h = 1.1221 

The prototype and model have different stacking sequence and number of the 
laminae (X N ± 1). Since the number of plies and their stacking sequence are not 
identical, then A^ is 


AtU CjX a A/j Xg 

where 

c = (512Q„+512Q 2a ) _ 3 

d (668.25Q„ + 175.5Q 22 ) * 

The experimental, predicted and theoretical values of the maximum deflection 
of the prototype and the experimental data of the model are presented in Figures 
4.7 and 4.8 . The predicted values of maximum deflection of the prototype are 
in excellent agreement with its experimental data. Figures 4.9 and 4.10 show 
the % discrepancy between the model and prototype. The % discrepancy between 
predicted and experimental results of the prototype are less than 5% while this value 
is 10% -15% for the theoretical and experimental data of the prototype, and 9% — 12% 
for theoretical and predicted data of the prototype. 
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Figure 4.9: %Discrepancy of Theory and Actual Test Results of Prototype G2 

( 0 / 90 / 0 /.. .)i 6 and its Model Gf (O3/9O3/O3/9O3/O3). 
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Figure 4.10: %Discrepancy for Prototype G2 ( 0 / 90 / 0 /...) is V Vhen Gf (O3/9O3/O3/9O3/O3) 
is Used as Model. 
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CASE - 3 


In this case, the E-Glass/Epoxy plate G4 is chosen as model which is used to 
predict the behavior of a Kevlar/Epoxy plate. The plate K7 is the prototype. 

From the data of Table 4.1 , the scale factors are calculated as; 

A 0 = 1.0 , A 6 = 0.9903 , A fc = 1.1832 

The prototype and its model have different numbers of plies with different stacking 
sequences. Similar to Case — 2, the A„, is, 


where 

(1584Q„ + 416Q 22 ) _ 3 

Cd (668.25g n + 175.5Q 22 ) N 

Figures 4.11 and 4.12 show the predicted and actual experimental deflections 
for different loads. The discrepancies for both the model and prototype are shown 
in Figures 4.13 and 4.14. The discrepancies between theory and experiment, and 
experiment and predicted for the prototype are very high (more than 35%). But the 
differences between theoretical and predicted data are less than 11%. Comparing 
the values of these two discrepancies ( experiment and theory, and experiment and 
predicted) of the prototype, they are in good agreement. 
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Figure 4.11: Predicted and Actual Test Results of Kevlar/Epoxy Plate K1 

(O4/9O4/O4/9CL1/O4) and its Model G4 (O3/9O3/O3/9O3/O3). 
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Figure 4.12: Theoretical, Predicted, and Actual Test Result of Prototype K 7 

(O4/9O4/O4/9O4/O4) By Using Model G\ (O3/9O3/O3/9O3/O3). 
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Figure 4.14: %Discrepancy for Prototype K7 (O4/9O4/O4/9O4/O4) When Model G4 

(O3/9O3/O3/9O3/O3) is Used. 
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CASE - 4 


In this case, the model and prototype have different stacking sequences of the 
laminae and ^ 1 . Plate K8 is the prototype and K6 is chosen as its scale model. 
From the data of Table 4.1 , the scale factors are calculated as; 

A 0 = 1.0 , A 6 = 0.9995 , A* = 1.0 

Since the number of plies and the stacking sequence of the laminates are not 
identical, then \ w is 


A nj c^A fl A* A, 


where 

(1584Q„+416g 22 ) _ 3 

Cd (729 Q n + 729 Q 22 ) N 

Similar to other cases, Figures 4.15 and 4.16 show the maximum deflections. 
The predicted results of the prototype match very well with the experimental data, 
especially for P > 600 lbs. Beyond this point the theory cannot predict the actual 
behavior of the prototype (possibly due to fiber or ply failure). Figures 4.17 and 4.18, 
which present the % discrepancies for the model and prototype, show this agreement 
between experimental and predicted data more clearly. As the load increases the 
discrepancy between theory and experiment for both model and prototype increase 
rapidly, while the discrepancies between predicted and experimental do not change 
significantly. 
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•’igure 4.15: Predicted and Actual Test Results of Kevlar/Epoxy Plate K8 

0 4 /90 4 /0 4 /90 4 /0 4 ) and its Model K6 (0/90/0/.. .)is 
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Figure 4.16: Theoretical, Predicted, and Actual Test Result of Prototype I\8 

(0 4 /90 4 /0 4 /90 4 /0 4 ) By Using Model K6 (0/90/0/. ..)i S . 
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Figure 4.17: %Discrepancy of Theory and Actual Test Results of Kevlar/Epoxy Plate 
K8 (O 4 / 9 O 4 /O 4 / 9 O 4 /O 4 ) and its Model K6 (0/90/0/. ..)i 8 . 
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Figure 4.18: %Discrepancy for Prototype K8 (O 4 / 9 O 4 /O 4 / 9 O 4 /O 4 ) When Model K6 

(0/90/0/. ,.) 18 is Used. 
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CASE - 5 


In this case, again the Kevlar/Epoxy plate K6 is chosen as the model and the 
Kevlar/Epoxy plate K9 is its prototype. From the data of Table 4.1 the scale factors 
are calculated as ; 

A. = 1.0 , X b = 0.9872 , A* =0.9184 

Since the prototype and its model have a different number of plies with different 
stacking sequences X w is 


Xyj 


Cd Xa 3 X h ~ 3 X g 


where 

_ (1404Qn j- 54 Q 22 ) . -3 

Q " (729 Q u + 729 Q 22 ) N 

Figures 4.19 and 4.20 show the predicted and actual experimental deflections for 
different loads. Similar to Case - 4 the predicted data of the prototype match very 
well with the experimental data, especially for P > 600 lbs. Beyond this point the 
theory cannot predict the actual behavior of the prototype (possibly due to fiber or 
ply failure). Figures 4.21 and 4.22 , which present the % discrepancies for the model 
and prototype, show this behavior clearly. As the load increases, the discrepancy 
between theory and experiment for both the model and prototype increase rapidly, 
while the discrepancy between predicted and experimental values does not change 
appreciably. The discrepancies for both model and prototype are shown in Figures 
4.21 and 4.22. 
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Figure 4.19: Predicted and Actual Test Results of Kevlar/Epoxy Plate K9 (06/906/0®) 
and its Model K6 (0/90/0/. ..)i 8 . 
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Figure 4.20: Theoretical, Predicted, and Actual Test Result of Prototype K9 

(0 6 /90 6 /0 6 ), By Using Model K6 (0/90/0/...)i*. 
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Figure 4.21: %Discrepancy of Theory and Actual Test Results of Prototype K9 

(0 6 /90 6 /0 6 ) and its Model K6 (0/90/0/.. .)is- 
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Figure 4.22: % Discrepancy for Prototype K9 (0 6 /906/0 6 ) When Model I\6 (0/90/0/...)x s 
is Used. 
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CASE - 6 


In this case the model and prototype have different stacking sequences of the 
laminae and A^r / 1 . Plate K9 is the prototype and K8 is chosen as its scale model. 
From the data of Table 4.1 the scale factors are calculated as; 


A a = LO , X b = 0.9877 , Xh = 0.9184 


Since the number of the plies and the stacking sequence of the laminates are not 
identical, then A^ is 


A^ — Cd A a A/j A, 


where 

_ (1404Q n + 54 Q 22 ) x —3 

(1584(5n + 4I6Q22) 

Similar to the other cases Figures 4.23 and 4.24 show the maximum deflections. 
The predicted values for the prototype match very well with the experimental data, 
especially for P < 600/6s. Beyond this point, the theory cannot predicted the actual 
behavior of the prototype. Figures 4.25 and 4.26, which represent the % discrepancy 
for model and prototype, show this behavior clearly. As the load increases, the dis- 
crepancy between theory and experiment for both the model and prototype increases 
, while the discrepancy between the predicted and experimental values increases at a 
slower rate. 
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Figure 4.25: %Discrepancy of Theory and Actual Test Results of Prototype K9 

(0 6 /90 6 /0 6 ) and its Model K8 (0 4 /90 4 /0 4 /90 4 /0 4 ). 
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Figure 4.26: %Discrepancy for Prototype K9 (0 6 /90 6 /0 6 ) When I<8 (0 4 /90 4 /0 4 /90 4 /0 4 ) 
is Used as Model. 


% Discrepancy of max. Deflection 


th. & pr. EM3 exp. 8c pr. 



100 200 300 400 500 600 700 


Load (lbs) 


54 





4.5 Stresses 


For the k th lamina, the stresses in terms of the strains and curvatures are; 

<5n Qn Q„ 1 '*' / r e°„ ) r k„ )\ 

= Q\2 Q 22 Q 24 [ < e°y> > + z < k yy m (4.29) 

. Q 14 Q 24 Q 44 J \ ( 1° xy J 1 k xy J / 

Where t° xx , e° yy , o.ndj° xy are the extensional and shear strains on the reference 
surface(z = 0) and k xx ,k yy ,k xy represent the change in curvature of the reference 
surface. The Qij are constant for a given lamina. In the case of cylindrical bending 
(u = u(x),v = 0, w = u>(x)) e 0 j, y = 7° xj/ = k yy — k xy = 0 which yields 




(4.30) 


(4.31) 


(4.32) 


°xx (k) = <?ff (u,* + ^W 2 X - ZW tXX 

Applying similitude theory for the normal stress, a xx , 

A...,, = **?<£ + £-*.£) 

'x /X x 


The resulting similarity conditions are; 


(4.33) 

(4.34) 


K xx w 


(4.35) 

K xx w 


(4.36) 

K xx w 


(4.37) 
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where A„, — A^AgA^j and A u — A^A^A^,, A^ 

To find which one of Eqs.( 4.35)— ( 4.37) gives the best prediction for the prototype 
behavior, the theoretical stress of the model is projected with each condition and 
compared to theoretical stress of the prototype. In this study, the theoretical stress 
(cr X x) of the model is considered as experimental stress of the model. This stress 
is projected by using derived similarity conditions, Eqs.( 4.35)-( 4.37), in order to 
predict the pertinent stress of the prototype. 

For complete similarity, Eqs.( 4.35) — ( 4.37) give the same result. However, for the 
distorted model each similarity condition gives different results. Figures 4.27— 4.34 
present the predicted and theoretical distributions of the normal stress in various 
layers of the prototype for cylindrical bending test. It is observed that the predicted 
stresses by Eq.( 4.37) agree very well with the theoretical results. Eq.( 4.36) cannot 
predict the behavior of the prototype accurately. Eq.( 4.35) is not a suitable similarity 
condition, since its predicted data do not match the theoretical results. The figures 
do not include the predicted stresses using Eq.( 4.35). The predicted stresses using 

Eq.( 4.36) are not included in all the figures. This is purposely done in order to 
simplify the figures. 
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Figure 4.27: Predicted and Theoretical Normal Stress a xx Distributions in Various 
Layers of the Prototype Gl (0/90/0...) 16 When G3 (O3/9O3/O3/9O3/O3) Is Used as 
Model. 
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Figure 4.28: Predicted and Theoretical Normal Stress a xx Distributions in Various 
Layers of the Prototype G2 ( 0 / 90 / 0 ...)i« Gf (O3/9O3/O3/9O3/O3) Is Used as 
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Figure 4 . 29 : Predicted and Theoretical Normal Stress cr xx Distributions in Various 
Layers of the Prototype K 7 (O4/9O4/O4/9O4/O4) When Gl ( 0 / 90 / 0...) 16 Is Used as 
Model. 



Figure 4 . 30 : Predicted and Theoretical Normal Stress cr xx Distributions in Various 
Layers of the Prototype K 7 (O4/9O4/O4/9O4/O4) When Gf (O3/9O3/O3/9O3/O3) Is Used 
as Model. 
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Figure 4.31: Predicted and Theoretical Normal Stress a xx Distributions in Various 
Layers of the Prototype K8 (O4/9O4/O4/9O4/O4) When K6 (0/90/0. ..)is Is Used as 
Model. 
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Figure 4.32: Predicted and Theoretical Normal Stress cr xx Distributions in Various 
Layers of the Prototype K9 (Oe/OOg/Oe) When K6 (0/90/0. ..)is Is Used as Model. 
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Figure 4.33: Predicted and Theoretical Normal Stress a xx Distributions in Various 
Layers of the Prototype K9 (0 6 /90 6 /0 6 ) When I<8 (O 4 / 9 O 4 /O 4 / 9 O 4 /O 4 ) Is Used as 


Figure 4.34: Predicted and Theoretical Normal Stress cr xx Distributions in Various 
Layers of the Prototype G5 (0/90/0...) 48 When Gl (0/90/0. ..) 16 Is Used as Model 
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